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METAL ARTIFACT CORRECTION IN COMPUTED TOMOGRAPHY 

Backgroun d ofthejfavention 

The following relates to diagnostic ima ^ ^ „ ^ ^ 
appltcahon ,„ computed tomography imaging of , ^ ^ incIudes hjgh 
such as meu, imptonts, dads, fiUh^ ». like> descrjbed ^ 

reference hereto. However , » also ^ ^ „ ^ ^ rf hfc 

(PET), three-dimensional x-ray imaging, and the like. 

to computed tomography imaging, soiled -metal artifacts" arise when the 
.magcd regton of interest contains mend implants, dental fiUings, or other regions of high 
radmnon absorption. Typ icall y, mem, artifacts appear in fte reacted image as ,JL 

.hanughm. region, of 

to* denatty are highly atftntmftd leading to substantia, measnremen, ereom, mtd the ffltered 
backprojection o, otiter reconahnction pmeess translates these measurement en-ora mm area* 
nnage artifacts commonly known as metal arti&cta. 

. . . „ A ^ onsm ^^»»««togme l al M Mfac te includes p«fom»ng fihered 

dt^TT te m mKomcM ■ »*• « 

densny tn ft. uncorrected reconrtnicftd image, and replacing projection fta. pass through Z 

vatoea. The coreecftd projection daft again undergoes filtered hackprojection to ptoduce a 
corrected reconstructed image. 

■..m ■ ^^^"^"^"'"^^'-^WficaHoninwhtchftere 
- a single, welWefined high density regta 8MTOm d e d „y much lower density tissue. It does 
network we,,, however, with a ptomli,, „f high density regions, or where there are medium 

me« are reduced but remain ve^ visible in fte cortecftd rented image 

especially be^een high density and medium deataity regiona. fc meto, imaging .ppficatioL 
medtum denarty regiona „ coneapond ft bmte while high 

cotreapond ft meta, imp,,*,, den,, ffllteg , opmtim «JLJU 
computed ftmography applications), „s devices, aatd the like. Hence, in ,J2 
~ ftmogmphy imaging, fte region of inftres. common* contain, memum^ 
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The present invention contemplates an improved apparatus and method that 
overcomes the aforementioned limitations and others. 

Brief Summary *> f the Invention 

According ,„ one aspect, a memod ia provided for producing , conMttd 
reacted tmage ftom acouired tomographic projection date. Acquired projection da., 
correspond** to a region are reconsSllcttd ^ „ ^ 

fte onco^d reconstruct image are ciassffleo a, ie« into high density, medmm density 
den,,,, and ,o„ density ctaasea m rep,aeed with pix e, ™,„es of the ■„„ density pix e, Cat 

IZL"£T T 110 ~ ta8e 15 *— * projecttd w • 

^"^•^^^ondateeonWbmmg,,,,.^,,^ 

n*Uced wuh corresponding aynthedc projeCon data to genem tt correct projecdon da«a. 

The conected section data are reconaaucKd into a co«cted reconstructed image 

According ,„ another aspect, an append is closed f „ produc ^ , 
coated record image from edited mmographic psojecuon data. A reconshJLg 
means . prodded for recmts.ruodng acquhed projecnon dara eormsponding » a region into an 
enc^reced recca^ image. A Cassaying me*t S „ prevMed for . £ ~ " 

ZTITT V m8e " ,eaM iMO U8h d ^ — d ~ - 

r IT. " enK °' mMnS " fM «*■*• Phtela of the tmcmreoteo 

for fo^rd projechng me aynftehc image „ genemte aynmeuc projecnon A projecnon 

P^*» data. The reconstmcdng meana reconsaucts me correct paction da,a m to a 
corrected reconstructed image. 

According to y„ another aspect, A mdiogmphic semmer ia diaoloeed including 
^■edWgr^ysc^^^^.^^^^^^ 

mount Tie scanner acquhee totnogmphic p^jeedon data during rotadoo ,f te „L g 

rsrzr -*~ a comoted image from *■ — «^ 

Zl d Perf0nni ° S ' " M ' h0d — *• acquired 

Proton d,,, eo^pc^dm, t „ a regi „„ mto an ^ « J 

■ow dens.,, p.,., Caasea; repiacing pixel, of me u.onecfcd reconshucted image ma, ^ of 
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the high density and low density classes with pixel values of the low density pixel class to 
generate a synthetic image; forward projecting the synthetic image to generate synthetic 
projection data- replacing acquired projection data contributing to the pixels of the high density 
class with corresponding synthetic projection data to generate corrected projection data; and 
reconstructing the corrected projection data into a corrected reconstructed image. 

One advantage resides in improved artifact correction for images having a 
plurality of discontinuous regions of metal or other high density material. 

Another advantage resides in improved artifact correction for images including 
a metal or other high density region and a bone or other medium density region. 

Numerous additional advantages and benefits will become apparent to those of 
ordmary skill in the art upon reading the following detailed description of the preferred 
embodiments. 



Brief Description « y the Drawings 
The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. The 
drawmgs are only for the purpose of illustrating preferred embodiments and are not to be 
construed as limiting the invention. 

FIGURE 1 diagrammatically shows a computed tomography imaging system 
including an artifact-correcting reconstruction processor. 

FIGURE 2 diagrammatically shows a preferred set of pixel density thresholds 
defining a preferred set of pixel density classes. 

FIGURE 3 shows a block diagram of the artifact-correcting reconstruction 
processor of the computed tomography imaging system of FIGURE 1. 

FIGURE 4 diagrammatically shows sinogram representations of acquired 
modified, and synthetic projection data at several points in the reconstruction processing. 

FIGURE 5 diagrammatically shows reconstructed and partially processed 
images at several points during the reconstruction processing. 

Detailed Descriptin,, of the Prefer red EmhnHim P „> c 
With reference to FIGURE 1, a helical conebeam computed tomography 
imaging scanner 10 includes a radiation source 12 that produces a radiation beam directed into 
an examination region 14. The radiation beam interacts with a region of interest of an imaging 
subject disposed in the examination region 14, producing spatially varying absorption of the 
radiation as it passes through the examination region. A radiation detector 16 detects the 
absorption-attenuated radiation after it passes through the examination region 14 
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In a preferred embodiment, the radiation source 12 produces a fan-beam or 
cone-beam of x-rays. The radiation source 12 and the detector 16 are preferably mounted in 
oppositely facing fashion on a rotating gantry 20 so that the detector continuously receives 
x-rays from the radiation source. As the source 12 and the detector 16 are revolved about the 
examination region 14 via the rotating gantry 20, views are acquired over an angular range of 
preferably about 360° or more. Optionally, a reduced scan of between about 180" and 360° is 
used, in another suitable embodiment, the detector 16 is replaced by a stationary detector ring 
mounted on a stationary gantry 22. Typically, a subject support 26 is linearly movable in an 
axial or z-direction that is generally transverse to a plane of rotation of the rotating gantry 20 

Multiple-slice computed tomography projection data are suitably acquired by 
performing successive axial scans with the subject support 26 stationary during each axial scan 
and stepped linearly between axial scan, In this arrangement, the detector 16 can have either a 
smgle row of detector elements (that is, a one-dimensional detector) or a two-dimensional array 
of detector elements. Alternatively, helical computed tomography projection data is suitably 
acquired during continuous linear movement of the subject support 26 and simultaneous 
rotatmg of the gantry 20. This effects helical orbiting of the radiation source 12 relative to an 
»magmg subject disposed on the subject support 26. A generally conical radiation beam and a 
two-Amensional radiation detector is preferably used for acquiring helical projection data 

The outputs of detector elements of the radiation detector 16 are converted to 
acquired integrated attenuation projection values ^ that are stored in an acquired projection 
data memory 30. Each projection datum ^ corresponds to a line integral of attenuation along 
a hne from the radiation source 12 to a corresponding one of the detector elements of the 
detector 16. The projection data can be represented in a sinogram format in which each 
^0-d.mensional slice of the imaged region of interest is represen|ed by a projection data array 
having coordinates of viewing angle (<j,) and line integral index (n). 

In the case of a fan-beam, cone-beam, or other original acquisition geometry 
havmg non-parallel rays, a rebinning processor 32 optionally rebins the projection data into 
parallel ^ews. For a cone-beam geometry, such parallel-rebinned views typically include some 
uncorrected beam divergence in the cone-angle direction. Although parallel rebinning can 
unprove computational efficiency, the artifact correction techniques described herein are 
typxcally most effective when performed in the original acquisition geometry, that is, without 
parallel rebmning. Hence, the parallel rebinning processor 32 is preferably omitted. 

For typical fan-beam and cone-beam geometries, the line integral index n 
suuably corresponds to a detector index indicating a detector element used to measure the 
project™ of index n. It is contemplated, however, that the line integral index n may lack a 
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direct corapondence w*h detector demen, oumber . Such a lack of direct com^dence ca n 
rcaulU»rcxa^ 1 c,tomm^oM OT b« W c„ rcbtae(fproJeMions . * 

aca«.l .r ^ ^ CyC ' taS P,0CeSS ° r 34 ^ tta,U8h 0,6 ""a™" ooncaponding to 

™* err' ^ -* bto - — — — — cL 

Pressor 40. Tie artrfact-co^g proMssor perfonm , 

denst* regto™ „ „„. or «. rocaiuca^ regfaB „ aMtio „ „ fc _ * 
density res™. The s^onoenng reconstruction _„ ^ rf . ™ £ 

^puce 1 c,as s ,amediu m aen sit yph t eIcU K , a , dalowdensil) , I)ixe , closs 

threads » ,IT T^T ^ ™ * " ' * « ° f *-» 

W « ■** ^ dasaes: , high denrfty ciass designated "A" contending to 

lass dented "B" cording » bone or other medfan, density features; a low denai* 
E corresponding » a. pockeu, ^ „, surromdfag „. J 

defiZ " <ta *" > ' " E "- A "** W-*m <taeity threshold H_„ 

dL^ ,OW 1 d r" yAra » S » ,M *-* *— H, define, a n^un, demity of * , w 
f° * ' "° ' "' aXi,Mra ° f *■ «— *» —» Caaa "D" A transit 

ctTJr T * fines a *-* ° f - — «- 

class D and a maximum density of the air density class "E", 

Although in the preferred embodiment five density classes "A", «B» - C « «D« 

mm the low dens.ty/trans.t.on density threshold H te is omitted and the threshold H a « , also 
«ie threshold 1S optlonal]y shifted fa a higher j 

transrtum p,xel class »D» between the low density pixel class "C" and th ■ a 

"F» fi,^u a • a «»"y pixel class C and the air density pixel class 

, S "" *" - »-> «*- » can ako he rnniBed, .eavtag J ^ ^ 

pixel classes "A" "B" »<-» *u ^ y aensity 

c . jI . A ' B * ° ' Wth *■ transitional and air density pixel classes »D« and »E» 
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-sa „ ,. Wi,h Mmi " e rcferenCe <° HGTJTtE I, for each tapu, slice the 
-option process 40 mrtpu* . tw^imensiou.! artilkotcoITCCled 
reconstruct Image. In m« ra . sltee or hefal 

zn^T^ ■"™ d — ~ — * - ^ ^ 

I" *T 8 ^i"™ «— — recited votame image, jf ^ 
ac q utred proj e«ion data „ limited „ „ ^ ^ rf fc ^ ^ 

ZlT*" *" c ™ tas 10 * 8 si " 8 ' 8 s,i ~ is *— - * - i 

™„o„ „ r <o and „e .mage ^ 46 ^ a ^^.J 
^<<— > recouped image. Optionally, p^tion ^ to om 

mme unage *«. are acquired over . ^ _ ta(mal (o provjde temporal seHes of 
artifcct-corrected reconstnicted imaiie slice? or i™„ . i 

„ , „ , . g OT mlage volumes representative of a temporal 

evolution of the region of interest. temporal 

memo™ « , " ^ Pr ° < * SSOr " """"^ "°~ ° r 3,1 ° f <* *. image 

«« ^ create a hmuan-vie^ble image representation sueh as a tae*^ 

lire, Tbe human-ve^bl. image representation is displayed on a user inters 52, Ln is 

^ dKPhym ' ^ * — « — * of image memory 4« car, * 

P™« on paper, stored in a non-volatile Cecbomc or magnetic smrag. medium transmit* . 

HTn ,! * «— *' computed tomography imaging 

se ss ,o»,n M mtorau,m,g,^ S e ss i < m, or „ theroiseopetatetllesoaiiiierlo ^ 

arofccw M nODRB 3 Sb<>WS 8 """* < "* en " n ° f « "^™« «nbodim«„, of me 
et^TT ™ C60n ~ 4 °- ™« " proton da, set P. «, 

correspondmg to a two^inmn^ slice of me region of interns, of me imaging sub act 
. Prefcraay, fc p[oj eenon data set « M udea potion da* ftom a 360- revoZn <^ 

™'l°l S dT " fte region 14: — — — - ~ 

^^«^«^ fe r^ 1 e,a b o« W angm. r cov^ i ,^ c ^ hted . 

me acouir d ^ "f 8 "" * * ' "preeenutn™ « of 

the acquired projection data set P. 60 is ikmm, ft.. . a- . 

„ ' " Sl ""™- Coordinate or y-eoordmate of the sinogram 

rep™,, « colTOponds to ^ ang]e fc fc ^ » 

smogram represen te do„ 62 correspond, to Hue m tegral ^ „. ^ ^ 

.2 t*,^*^^ re8io „ ^z^ZZ 
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of a medium density region. In general, the slice can include an arbitrary number of high 
dens.ty regions and an arbitrary number of medium density regions, as -well regions of other 
classes shown in FIGURE 2. The regions produces various corresponding sinogram trajectories 
that can cross one another. 

With returning reference to FIGURE 3, a filtered backprojection processor 70 
performs filtered backprojection of the projection data set P„ 60 to generate an uncorrected 
reconstructed image I 74. Although filtered backprojection is preferred, the processor 70 can 
alternatively implement substantially any type of image reconstruction algorithm that is 
compatible with the geometry of the projection data set P 0 60. 

With reference to FIGURE 5, as is known in the art, the presence of one or 
more mgh density regions in the slice typically causes the uncorrected reconstructed image I, 
74 to include metal artifacts which generally manifest in the image as streaks extending away 
from the high density region. FIGURE 5 shows a diagrammatic representation of an exemplary 
uncorrected reconstructed image I. 74 which diagrammatically shows such streaks The 
d,agrammatic reconstructed images of FIGURE 5 are not intended to correspond to images of 
any particular anatomical structure, but rather diagrammatically illustrate features of the 
arhfact-correcting reconstruction process that are typically observable in image space. 

With returning reference to FIGURE 3, pixels of the uncorrected reconstructed 
.mage 1. 74 are classified by a pixel density classification processor 78 to generate a segmented 
or p,xel density-classified image W 80 in which pixel values are replaced by density 
classification index values corresponding to the high, medium, low, transition, and air density 
Pixel classes of FIGURE 2. In a preferred embodiment, the pixel density classification 
processor 78 employs thresholding using the set of thresholds [H^ 42 to classify each pixel of 
the uncorrected reconstructed image I, 74 into an appropriate pixel density class. 

The pixel density-classified image W 80 is diagrammatically shown in 
FIGURE 5, where regions of the image essentially consisting of pixels of a particular class are 
labeled by an appropriate class index selected from the density class indices "A" «B» »C» 
»D" "E" of FIGURE 2. As indicated in FIGURE 5, me exemplary image includes a region of 
mgh density »A» which could be a metal implant and a region of medium density »B« which 
could be a region of bone, both contained within a region of low density »C» which could be 
soft t,ssue or the like. An annular region of transition density «D« surrounds the region of low 
density »C», and a region of air density »E» fills the periphery of the pixel density-classified 
•mage U. 80, corresponding for example to ambient air surrounding the imaging subject. 

Although the exemplary transition region in FIGURE 5 is an annular boundary 
regmn, the transition density class »D« is not limited to annular boundary regions, the 
thresholds H Ic and H MttBX) defining the transition density class «D» are preferably selected so 
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that the transition density class »D« includes air cavities and regions surrounding air cavities 
inside human subjects. For example, the transition density class «D« preferably includes air 
. cavities in the head such as sinus cavities, the throat, nasal cavities, and the like. 

With returning reference to FIGURE 3, a density class averaging processor 84 
computes an average density value for pixels of the low density class »C». m preferred 
embodiments which include an air density class, the density class averaging processor 84 also 
optionally computes an average density value for pixels of the air density class »E» Averaging 
of the air density class is omitted in the embodiment that employs only three classes and omits 
the air density class altogether. 

A correction image processor 88 selectively replaces selected higher density 
P«ek of the uncorrected reconstructed image J. 74 with lower density values to produce a 
synthetic image W 90. Specifically, pixels of the uncorrected reconstructed image 1. 74 
ailing within the high density class "A" as indicated by the pixel density-classified image W 
80 are replaced by the average density value for pixels of the low density class »C». Similarly 
pixels felling within me low density class „ c „ ^ fey ^ ^ ^ ^ 

Pixels of the low density class »C». Alternatively, the »A« and »C» class pixels can be replaced 
with a preselected value or spectrum of values which fall near the center of a nominal »C» 
density region. If an average density value for pixels of the air density class has been 
computed, pixels falling within the air density class "E" are suitably replaced by the average 
density value for pixels of the air density class "E". Pixels of the medium density class »B» are 
not replaced. If a transition density class »D» has been defined, then pixels of the transition 
density class "D« are also preferably not replaced. 

With reference to FIGURE 5, the synthetic image 1^ 90 is diagrammatically 
shown. It will be observed that the synthetic image W 90 contains contrast principally due to 
the regions of medium density -B- and of transition density "D". Regions of high density "A" 
are effectively removed from the synthetic image W 90 by replacing pixels in these regions 
with the average density value for pixels of the low density class »C». Similarly, image contrast 
due to density variations in regions of low density »C» are substantially reduced by substitution 
of these pixels by the average density value for pixels of the low density class ''C» 

FiriTPP * , retUmin8 10 nGURE 3 and ^ additional refe ™«* to 

FIGURE 4, a forward projection processor 94 forward projects the synthetic image W 90 to 

produce synthetic projection data Psymh 96. A diagrammatic sinogram representation 98 of the 

projection data set 96 is shown in FIGURE 4. The sinogram representation 98 retains the 

exemplary trajectory 66 of the medium density region of the sinogram representation 62, since 

the Pixels of medium density Cass «B" were not replaced. Similarly, trajectories due to regions 

of transrtion class »D« are retained in the synthetic projection data P^ 96 since these pixels 
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are retained in the synthetic image ^ 90. The sinogram representation 98 does not include 
the exemplary trajectory 64 of the high density region of the sinogram representation 62 
however, since the pixels of high density class "A" were replaced by the average density value 
for pixels of the low density class "C". 

In general, the sinogram representation 98 of the synthetic projection data set 
P sywb 96 retains attenuation contributions to the line integrals due to regions of medium density 
class »B«, even if those line integrals additionally have contributions from regions of high 
density class «A» in the acquired projection data set P 0 60. The forward projection processor 94 
preferably projects the synthetic image 1^ 90 using a projection geometry corresponding to 
the geometry of the computed tomography imaging scanner 10. 

With continuing reference to FIGURE 3, the synthetic projection data set P^ 
96 .s used for selectively replacing projection data of the acquired projection data set P 0 60 
having attenuation contributions from high density regions. A high density region edge finding 
processor 100 identifies edges of regions essentially consisting of pixels of the high density 
Pixel class. Regions of high density can be identified, for example, by employing a moving 
analysis window and identifying high density regions as those having more than a selected 
number of pixels of the high density class within the analysis window. Once a high density 
region ,s identified, a suitable edge finding algorithm employs image differentiation or 
convolution to selectively enhance and identify edges. 

In an alternative approach, the edge finding processor 100 performs binary 
thresholding on the uncorrected reconstructed image I, 74 using ^ high density/medium 
density threshold to produce a binary image having "1" binary values for pixels of the 
high density class and "0" binary values for other pixels. The binary 'T values are then filtered 
to remove outlying pixels of value »1« that have few or no neighboring pixels of value »1» and 
the remaining groups of pixels having value »1« define the high density regions. Edges are then 
identified in the filtered binary image as »0«-to-«l« binary value transitions, identified for 
example using an "exclusive-or" binary operator acting on neighboring pixel pairs. Those 
sklled in the art can readily employ other edge finding algorithms in constructing the edge 
finding processor 100. 

The edge finding processor 100 suitably outputs an edge-enhanced image 
Iho.^ 102 that identifies edges of high density regions. The edge-enhanced image 1^ ^ 102 
is diagrammatically shown in FIGURE 5, where edges 106 of the region of the uncorrected 
reconstructed image I. 74 consisting essentially of pixels of the high density class are indicated 
In exemplary FIGURE 5, the edges 106 define a substantially circular edge; however, the high 
density regions can have substantially arbitrary shape with substantially arbitrary edges 
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Typically, the high density regions correspond to dental fillings, metal implants, or other 
compact discrete objects, and have well-defined edges defining simply-closed geometries. 

With reference to FIGURE 3, a high density pixel replacement processor 110 
selectively replaces projections of the acquired projection data P 0 60 having contributions from 
high density regions with synthetic projection data from the synthetic projection data set P syitth 
96 to produce a corrected projection data set P, 112. In one approach, the high density pixel 
replacement processor 110 replaces projection data that intersect one or more high density 
regions identified by the edge-enhanced image V.^ 102 with corresponding projection data 
from the synthetic projection data set P^ 96. This approach, however, can lead to substantial 
attenuation discontinuities at transitions between the original acquired projection data P„ 60 and 
the replacement synthetic projection data set P sytlth 96. 

Hence, in a preferred embodiment, the high density pixel replacement 
processor 110 interpolates between the acquired projection data and the synthetic projection 
data at the identified edges 106 of the high density image region. A suitable interpolating 
projection replacement formula replaces the acquired projection data of the projection data set 
P„ 60 with replacement projection data ud^, having values given by: 1 



^ repl (n)=^d synm (n) + af-5^ > | + bf n ~ n ' 1 



(1), 



where a=[ndo(n,)- U<Wn,)], b=[nd„(n 2 )- ud syMh (n 2 )], index n is the line integral index, indices 
n, and n 2 are line integral indices of the edges 106 of the high density image region as shown in 
FIGURE 5, tido indicates acquired projection data of the projection data set P„ 60, and ud 8ynUl 
indicates projection data of the synthetic projection data set P^ 96. Projection data in the 
range n,<n<n 2 is replaced in accordance with Equation (1). Review of Equation (1) shows that 
^i(n,y= udo(n,) and ud^fa^ udo(n 2 ), providing a smooth transition at the edges 106. 

The interpolative replacement performed by the high density pixel replacement 
processor 110 operates on each view of the sinogram that includes one or more high density 
regions. That is, for each view specified by a given viewing angle <>, a given high density 
region has edges at line integral indices n„ n 2 which in general are different for different views 
of the same high density region. The interpolative replacement is performed for each view 
using the n„ n 2 edge values computed for that view. Moreover, there may be more than one 
non-contiguous high density region within a given view, with each such non-contiguous high 
density region having its own n„ n 2 edge index values. The interpolative replacement set forth 
in Equation (1) is repeatedfor each non-contiguous or non-overlapping high density region. 
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With reference to FIGURE 4, a diagrammatic sinogram representation 114 of 
the corrected projection data set P, 112 is shown. Advantageously, the sinogram representation 
114 retains the trajectory 66 of the medium density region substantially intact, even where the 
trajectory 66 of the medium density region crosses the location of the removed trajectory 64 of 
a high density region (indicated in the diagrammatic sinogram representation 114 by dashed 
hnes representing the edges n„ n 2 of the replaced high density region in the views). There are 
no gaps of low density projections along the trajectory 66 of the medium density region 
because the synthetic projection data set J> syn<b 96 retains contributions of medium density 
regions to the integrated attenuation projection values ud„ while selectively removing 
contributions of high density regions to the integrated attenuation projection values ud, 
Similarly, although not illustrated, a trajectory of a transition density region «D« remains intact 
even where such a trajectory is crossed by the high density trajectory 66. 

The corrected projection data set P, 112 is input to the filtered backprojection 
processor 70 which performs filtered backprojection to generate an artifact-corrected 
reconstructed image 120. Preferably, a labeling processor 122 substitutes pixels defining a 
preselected label for pixels of the artifact-corrected reconstructed image 120 corresponding to 
h.gh densny regions that were substantially modified by the substitution of interpolated 
synthetic data. This labeling notifies a radiologist or other user who views a final 
artrfact-corrected reconstructed image 124 that the labeled image regions are substantially 
modrfied by the artifact-correcting process. The pixels defining the preselected label can for 
example, be pixels having high density values simulating the original high density region. In 
another approach, the pixels defining the preselected label define a selected pattern, such as a 
cross-hatch pattern, that is clearly artificial. The final artifact-corrected reconstructed image 
124 is stored in the image memory 46 of FIGURE 1 . 

The invention has been described with reference to the preferred embodiments 
Obviously, modifications and alterations will occur to others upon reading and understanding 
the preceding detailed description. It is intended that the invention be construed as including 
all such modifications and alterations insofar as they come within the scope of the appended 
claims or the equivalents thereof. 
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Claims 



Having described the prefeired embodiments, the invention is now claimed to be: 

1. A method for producing a corrected reconstructed image from acquired 
tomographic projection data, the method including: ' 

reconstructing acquired projection data (60) corresponding to a region into an 
uncorrected reconstructed image (74); 

classifying pixels of the uncorrected reconstructed image (74) at least into high density, 
medium density, and low density pixel classes; 

replacing pixels of the uncorrected reconstructed image (74) that are of the high density 
and low density classes with pixel values of the low density pixel class to generate a synthetic 
image (90); 

forward projecting the synthetic image (90) to generate synthetic projection data (96); 

replacing acquired projection data (60) contributing to the pixels of the high density 
class with corresponding synthetic projection data (96) to generate corrected projection data 
(112); and ( 

reconstructing the corrected projection data (112) into a corrected reconstructed image 

(120). 

2. The method as set forth in claim 1, wherein the replacing of pixels of the 
uncorrected reconstructed image (74) that are of the high density and low density classes with 
pixel values of the low density pixel class includes: 

computing an average value of pixels of die low density pixel class; and 
replacing pixels of the uncorrected reconstructed image (74) that are of the high density 
and low density classes with the average value of pixels of the low density pixel class. 

3. The method as set forth in claim 1, wherein the region is a slice and the uncorrected 
reconstructed image (74) is a two-dimensional image, the method further including: 

repeating the reconstructing of acquired projection data (60), the classifying, the 
replacing of pixels, the forward projecting, the replacing of acquired projection data (60) and 
the reconstructing of the corrected projection data (112) for acquired projection data 
corresponding to each of a plurality of two-dimensional slices; and 

combining the two-dimensional reconstructed images (120) corresponding to the 
plurahty of two-dimensional slices to generate a three-dimensional corrected reconstructed 
image. 
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4. The method as set forth in claim 1, wherein the replacing of acquired projection 
data (60) contributing to the pixels of the high density class with corresponding synthetic 
projection data (96) includes: 

during the replacing, interpolatively adjusting the synthetic projection data (96) to 
smooth transitions between the synthetic projection data (96) and the acquired projection data 



(60). 



data J: ™T eth0d " fortH " Cl3im Wh<!rein rCplaCin8 ° f aC ^ ed Potion 
data (60) contrabutmg to the pixels of the high density Cass with corresponding synthetic 

projection data (96) includes: 

identifying a high density image region substantially comprised of pixels of the high 
density class; and 

performing the replacing for acquired projection data (60) contributing to the high 
density image region. 

6. The method as set forth in claim 5, wherein the region is a slice, the uncorrected 
reconstructed image (74) is a twoKlimensional image, the acquired projection data (60) is in a 
smo^am^ 

the hrgh densrty class with corresponding synthetic projection data (96) further includes- 

replacmg the acquired projection data (60) with replacement projection data pdL, 
havmg values given by: M ^ repl 

(n) = Md^ (n) + af -HlZlO + b f p ~°. ) 

where a^n,)- ud^n,)], b=[ud <) (n 2 )- ud^], n is a line integral index, indices n, and 
n 2 are hne integral indices corresponding to the edges (106) of the high density image region 
with replacement being performed for n,<n<n 2 , ^ indicates acquired projection data (60), and 
JWsynfc indicates synthetic projection data (96). 

7. The method as set forth in claim 5, wherein the region is a slice, the uncorrected 
reconstructed image (74) is a two-dimensional image, the acquired projection data (60) is in a 

the mgh densuy class with corresponding synthetic projection data (96) further includes: 
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for at least one view of the sinogram corresponding to a view angle, repeating the 
KJenUfymg of a high density image region and the performing of the replacing for acquired 
projection data (60) contributing to the high density image region for at least two 
non-contiguous high density image regions. 

8 The method as set forth in claim 5, wherein the identifying of a high density image 
region substantially comprised of pixels of the high density class includes- 

applying an edge detection algorithm to identify edges (106) of the high density image 
region. & 



9. The method as set forth in claim 8, wherein the performing of the replacing for 
acquired projection data (60) contributing to the high density image region includes- 

datarJT^ bCtWeen ^ a ° qUired Pr ° jeCti0n ^ (60) *• 
data (96) adjacent the identified edges (106) of the high density image region. 

10. The method as set forth in claim 5, wherein the identifying of a high density image 
region substantially comprised of pixels of the high density class includes: 

. fc ff™** bina <* ^holding of the uncorrected reconstructed image (74) using a 
threshold that assigns a first binary value to pixels of the high density pixel class and that 
assigns a second binary value to pixels not of the high density pixel class; and 

identifying edges (106) of the high density image region as transitions from pixels of 
the first binary value to pixels of the second binary value. 

11. The method as set forth in claim S, further including: 
superimposing a label on the high density image region. 

12 The method as set forth in claim 1, wherein the reconstructing of the acquired 
projection data (60) and the reconstructing of the corrected projection data (112) each include- 
reconstructing the projection data (60, 112) using filtered backprojection. 

data ,60V fT meth ° d 38 f0rth " ^ Wherdn * e aCqUirCd •—W* P*— 
data (60) is from a computed tomography scanner (10) and is in a cone-beam geometry, and the 
forward projecting includes: 

forward projecting the synthetic image (90) to generate synthetic projection data (96) 
in the cone-beam geometry. 
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14. The method as set forth in claim 1, wherein: 

the high density pixel class corresponds at least to metallic material; 
the medium density pixel class corresponds at least to bone; and 

the low density pixel class corresponds at least to soft tissue 

l 

15. The method as set forth in claim 14, wherein: 

the classifying of pixels of the uncorrected reconstructed image (74) further includes 
classifying pixels of the uncorrected reconstructed image (74) into an air density pixel class 
having a maximum density that is lower than a minimum density of the low density pixel class; 

the replacing of pixels of the uncorrected reconstructed image (74) to generate the 
synthetic image (90) includes replacing pixels of the uncorrected reconstructed image (74) that 
are of the air density class with an average value of pixels of the air density pixel class. 

16. The method as set forth in claim IS, wherein the classifying of pixels of the 
uncorrected reconstructed image (74) further includes classifying pixels of the uncorrected 
reconstructed image (74) into a transition density pixel class spanning a density range between 
the maximum density of the air density pixel class and the minimum density of the low density 
pixel class. J 

17. The method as set forth in claim 1, wherein the replacing of pixels of the 
uncorrected reconstructed image (74) does not include replacing pixels of the medium density 
pixel class. ' 

18. An apparatus for producing a corrected reconstructed image (120, 124) from 
acquired tomographic projection data (60), the apparatus including: 

a reconstructing means (70) for reconstructing acquired projection data (60) 
corresponding to a region into an uncorrected reconstructed image (74); 

a classifying means (78) for classifying pixels of the uncorrected reconstructed image 
(74) at least into high density, medium density, and low density pixel classes; 

a pixel replacement means (88) for replacing pixels of the uncorrected reconstructed 
-mage (74) that are of the high density and low density classes with pixel values of the low 
density pixel class to generate a synthetic image (90); 

a forward projecting means (94) for forward projecting the synthetic image (90) to 
generate synthetic projection data (96); and 
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a projection replacement means (100, 110) for replacing acquired projection data (60) 
contnbutmg to the pixels ofthe high density class wim corresponding acetic projection data 
1*0) to generate corrected projection data (112)- 

M T ft,di ° s (I0,) for **■ < 1M) of - - — «* 

I e r onsco^s ta8esK „,i, l l y „fp ixels(>fthehighdensitj , ])ijtelcUss * 

(100, 110^'form^uichKias: 38 ^ ' B c ' a ' nl w herein the projection replacement means 

an interpolating replacement means (H0) for interpolatively replacing acquired 
~ data „ with ^ ^ ^ ^ a , fte ^ ( ^ J fc 

21. The apparalua as set forth in claim 19, tether molndtag- 
regfol recons^ted image (120, corresponding fo the on, or .ore high density 

a set of T HZ "T"" " " ** " ^ 18 ' " herei " *" "-a** — . (78) accesses 
a set of threshold values (42) including at least: 

a medten tfensiry/high density threshold value defining a minimum 
dens,* „f fte Ush denaaypise, class md , amhama 
density pixel class, and 

a low density/medium density threshold value defining a rninimum 
densny ofthe medium density pixel class and a maximum density ofthe low 
density pixel class; and 

l*^" C,8SS ^ *"* <* •» — « mconaoncted image (74) a, 
fc- mfo htgh density, medta, ^ , ow ^™ « 

pixel values with the set of threshold values (42). "mpartson of 

tefcerin" ~ " ^ fcrthfa " - va-ues (42, 
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a low density/transition density threshold defining a minimum density 
of the low density pixel class and a maxi mum densily of a ^ densi 
pixel class, and 

a transition density/air density threshold defining a minimum density 
of the transition density pixel class and a maximum density of an air density 
pixel class, 

-» (* otto transition denshy ^ * densjty pixa J£ 
values with the set of threshold values. P 

« The apparatus as set for* to chdtn 23, „ hereto „. pixel 

M T 6 '° ,CO,reC,ed ""^ ™ *- dely 

pncel class rvtfo an avenge value of pbtels of foe air density pixel elass. 

25. The apparatus as set forth fa olaim 18. further including- 

*™-***«««--«*--.«> .adetector array < 16 >, a*, a rotafo,g gantry 

toll? ~" C,2> ^ ^ <16 > - —4 -he ZZ 

tomegtaphtc proton data (60, being ^ ^ „ fc ^ ^ ^ 

26. A radiographic scanner including: 
array*^^ 

quired ZTn ™ ^ Pr0dUCiD8 3 " ™ CtCd ta "* < 120 ' 124 > *» 

reconstructing acquired projection data (60) corresponding to a region 
into an uncorrected reconstructed image (74), 

classifying pixels of the uncorrected reconstructed image (74) at least 
into high density, medium density, and low density pixel classes, 

replacingpixelsofthe uncorrected reconstructed image p4) that are of 
the h,gh density and low density classes with pixel values of the low density 
pixel class to generate a synthetic image (90), 

forward projecting the synthetic image (90) to generate synthetic 
projection data (96), 
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replacing acquired projection data (60) contributing to the pixels of the 
high density class with corresponding synthetic projection data (96) to generate 
corrected projection data (1 1 2), and 

reconstructing the corrected projection data (112) into a corrected 
reconstructed image (120). 
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METAL ARTIFACT CORRECTION IN COMPUTED TOMOGRAPHY 



/ 

Abstract 



An artifact-correcting image reconstruction apparatus includes a reconstruction 
processor (70) that reconstructs acquired projection data (60) into an uncorrected reconstructed 
»mage (74). A classifying processor (78) classifies pixels of the uncorrected reconstructed 
-age (74) at least into high, medium, and low density pixel classes. A pixel replacement 
processor (88) replaces pixels of the uncorrected reconstructed image (74) that are of the high 
density and low density classes with pixel values of the low density pixel class to generate a 
synthetic image (90). A forward projecting processor (94) forward projects the synthetic image 
(90) to generate synthetic projection data (96). A projection replacement processor (100, 110) 
replaces acquired projection data (60) contributing to the pixels of the high density class with 
corresponding synthetic projection data (96) to generate corrected projection data (112) The 
reconstruction processor (70) reconstructs the corrected projection data (112) into a corrected 
reconstructed image (120). 
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